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ABSTRACT 

Extremely relativistic particles play an important role in various 
radiative mechanisms of pulsar radio as well as high-energy emission. It 
is thus very essential to estimate the Lorentz factor of these particles 
in the radiation regions both observationally and theoretically. In this 
paper, according to the frequency dependence of component separation 
and considering some radiation processes, we obtain the Lorentz factor 
and its variation in the radiation region observationally for seven 
multi- frequency observed pulsars. It is found that the Lorentz factor 
decays substantially in the emission region, and the decay can not be 
well understood by all known energy loss processes. This result hints 
that there are some unknown processes which are responsible to loosing 
energy effectively. 

Subject headings: Pulsars - Radiation mechanisms 
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1. Introduction 

The energy loss of relativistic particles in pulsar magnetosphere have been 
studied by many authors (Daugherty & Harding 1989, Sturner 1995, Zhang et 
al. 1997, Lyubarskii & Petrova 2000, hereafter LP2000). Three possible energy 
loss mechanisms, i.e., inverse Compton scattering (ICS) of thermal photons by 
relativistic particles, triplet pair production and curvature radiation of particles, 
were investigated by Sturner (1995). It is found that the resonant ICS process in 
strong magnetic fields is the dominant energy loss mechanism when the Lorentz 
factors are between 10 and 10 4 and pulsar surface temperature is ~ 10 6 K. 
Considering resonant ICS of thermal photons by secondary particles (10 < 7 < 10 4 ) 
above the pulsar inner gap, LP2000 found that the ratio of Lorentz factor to its 
initial value, 7/70, decreases generally to 40% within a height <; 10 km along 
the magnetic field lines above the polar cap and nearly keeps constant at higher 
distances. 

The Lorentz factor 7 can be estimated from some observational properties (e.g. 
microstructure width) of pulsars. Lange et al. (1998) assumed that the widths of 
microstructure are related to the beam widths of relativistic particles. They found 
that the derived 7, according to the observed micropulse typical widths of PSR 
B0950+08, PSR Bl 133+16 and PSR B2016+28, is in a range from 240 to 520 if 
the angle between the line of sight and pulsar rotation axis is 90°. 7 could be even 
larger since this angle may be less than 90° (Lyne & Manchester 1988, hereafter 
LM88, Rankin 1993). 

However, as far as we know, there are few works trying to derive energy 
loss of particles in the emission region observationally. Thorsett (1991, hereafter 
T91) studied the dependence of component separation on observing frequency 
for seven pulsars. It is found that the separation of emission components varies 
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with observing frequencies smoothly. According to this observational analysis and 
employing pulsar radio emission models, we calculate the 7 factors at different 
heights. It is found that in the radio emission region, the Lorentz factor decays 
significantly, which can not be understood by the energy loss processes known 
hitherto. 

In §2 the observational facts of the frequency dependence of component 
separation are introduced. The basic assumption, calculation and results are 
showed in §3. In §4, conclusions and discussions are presented. 

2. The frequency dependence of component separation 

It is based on the dependence of component separation AO on observing 
frequencies v to find out the energy distribution of particles in the emission region 
of pulsars. The dependence has been studied for a long time. Earlier analysis found 
that a simple power law could not fit the frequency dependence of component 
separation (hereafter FDCS) for most pulsars (Lyne et al. 1971, Sieber et al.1975, 
Rankin 1983, Slee et al. 1987). Traditionally two power laws were introduced, 



where the break frequency is usually between 0.14-1.5GHz, a<i < a\ indicating 
that at higher frequencies the component separation is less frequency dependent. 
Rankin (1983) even introduced a third power-law at intermediate frequency to fit 
PSR B1133+16. Collecting all available data, T91 restudied the FDCS for seven 
pulsars that exhibit double- or multiple-component average profiles, and proposed 
a smooth function, 




(1) 



AO = Av~ n + A0 mi 



(2) 
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to model for each pulsar (for multiple-component profiles, the component separation 
is the width between the two peaks of outmost cone). As T91 pointed out, it was 
mainly the insufficient sampling that led those earlier works to report a break 
frequency. Similar to the FDCS, the average profile width measured at 50% or 10% 
of component peaks also narrows as frequency increases (Rankin 1983, Xilouris et 
al. 1996). 

Our following analysis and results are based on Eq.(2). On the one hand, the 
index n and minimum component separation A# m j n were given by T91. On the 
other hand, Ad is a geometrical result which depends on emission height, and v is 
model-dependently related to the Lorentz factor in the emission region. Combining 
these two hands, we can derive a relationship between the distances and the Lorentz 
factors, when the geometric effect and a radio emission model are involved. 



3. Results: Electron Lorentz factor and emission distance 

Two assumptions are used in our analysis: 1) the magnetic field lines are 
dominantly dipolar in pulsar radio emission region, 2) radio emission arises from 
the last open field lines. The beaming angle 9^ (the angle between the tangent 
of magnetic field and the magnetic axis) can then be calculated from the viewing 
geometry (LM88), 

sin 2 (y ) = sin2 sin a sin ( a + + sin 2 (^j , (3) 

where the impact angle £ is the angle between the line of sight and magnetic axis, a 
is the inclination angle, a and C can be obtained from observations (LM88, Rankin 
1993). Geometrically, one has the following equations (Qiao & Lin 1998, hereafter 
QL98) about 9^, from which the distance r can be calculated, 

3 sin 9 cos 9 

tan ^=(2-^)' (4) 



-6 - 



sm 2 9 = —, (5) 

where the polar angle 9 is the angle between position vector r (from neutron star 
center to emission point) and the magnetic axis, r = |r|, i?LC = is the radius of 
light cylinder, with c the light velocity and P the pulsar rotating period. 

The Lorentz factor 7 of relativistic particles can be model-dependently 
obtained for a given r. Among pulsar radio emission models, the curvature 
radiation model (e.g., Ruderman & Sutherland 1975, hereafter RS75) and the ICS 
model (Qiao 1988, QL98) predict different functions of v = u(-y,r), which will be 
used below to calculate 7 = 7(7*), respectively. 

In the RS75 model, relativistic particles give out emission through curvature 
radiation with frequency 

- = h 3 CR ( , (6) 



2 \2-7rp 
then, 

7CR=(^ r j , (7) 

where 'CR' denotes 'curvature radiation', p is the curvature radius at a given point, 
which reads 

_ (l + 3cos 2 0) 3 / 2 sing 
P ~ 3(l+cos 2 #)# LC " () 

In the ICS model, low frequency electromagnetic wave, with frequency of 
uq ~ 10 6 Hz, is produced by inner gap sparking (though the frequency may vary 
for different pulsars due to different gap parameters, it does not change the results 
essentially for a given pulsar, so we simply take the frequency as v$ ~ 10 6 Hz). 
These low-frequency photons are up-scattered by relativistic electrons to higher 
frequencies (QL98), 

3 

v = -'tf GS v Q (\-Pca&Oi), (9) 
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then one has 

7ICS = UJ (10) 
where /? = (1 — I/t^cs) 1 ^ 2 ) the incident angle B\ is the angle between the wave 
vector of low frequency wave and the direction of electron moving along field lines, 
which can be calculated to be (QL98), 

CQze _ 2cos6>+ {R/r)(l -3cos 2 fl) 

COS 1 ~ y/(l + 3 cos 2 W) [1 - 2{R/r) cos 6 + {R/rJ 2 \ ' 1 j 

Employing the observational parameters listed in column 2 to 7 in Table 1, 
viz., the frequency range, A, n, and A6 m - m given by T91, a and £ given by LM88, 
we calculate 7 and r numerically and present the data in Fig. 1 for each pulsar. 
The values of minimum distance r m i n , maximum Lorentz factors ^ycR,max 

and 

7CRmax, and the ratio of 7CR " min and 7ics ' min are listed in table 1, where r min 

' ^min 7CR,max 7lCS,max 

and 7cR,max(7ics,max) are obtained from the maximum observable frequency, while 
r max and 7cR,min(7ics,min) are from the minimum observable frequency. 



4. Conclusions and discussions 

According to the frequency dependence of component separation of seven 
pulsars (T91), we obtain a model-dependent relationship between the Lorentz 
factors 7 of extremely relativistic particles and the emission distances r for each 
pulsar. From the analysis, we have following two points. 

1) . The radio emission of outmost cone typically arise from the distance of 
about 100km to about 300km, which is consistent with the previous results for 
these pulsars (e.g. Phillips 1992, Hoensbroech & Xilouris 1997). 

2) . Either 7cr or 710s decreases monotonically as r increases. Within the 
range of r, 7cr may reduce by about 60% and 90 % with respect to their maximum 
values, respectively. The results reflect that particles undergo severe energy loss 



-8 - 



in the radio emission region, which is essentially inconsistent with the theoretical 
prediction made by many authors, for example, LP2000, as will be discussed below. 

Some discussions related to our analysis are as follows. 

1) . In our analysis, emission is produced at the last open field lines. In 
fact, to derive 7(r), one can also assume that radio emission arises from inner 
magnetic fields, and the results are nearly the same, for an example, with the 
radius R c = XRlc, where the factor A > 1 (this means that the field lines are closer 
to to the magnetic axis, see QL98). Fig. 2 compares the curves of 7(7") obtained 
for the last open field lines and inner field lines with A = 1.2 (corresponds to the 
field lines satisfying 6 = 0.916>p, where dp is the polar cap angle defined by the last 
open filed line at pulsar surface) for PSR B 0525+21, from which one can see that 
the j(r) does not show significant difference. The uncertainty of inclination angle 
a and impact angle £ may influence the 7(7*), as shown in Fig. 3, which presents 
the 7(7*) curves for three groups of a and £ , i.e., (23.2°, 0.7°) (given by LM88), 
(23.2°, 2.0°) and (90°, 0.7°). Here It shows that the uncertainty of a may change 
the results of r greatly, but in general, 7(7-) does not change essentially. 

2) . Interestingly, our results show that there must be some additional 
energy-loss processes in the radio emission regions. But so far we did not find any 
reasonable effective mechanism for the energy loss. Theoretical study, e.g., LP2000, 
shows that 7 nearly does not decay at the emission region higher than 10km if only 
electron ICS process of the thermal photons from neutron star surface is included. 
Also pulsar radio emission can not lead electrons loose energy significantly (Lin 
1997). 

A possible effective process may be the resonant ICS of soft X-rays produced 
in outer gap with the relativistic particles. Cheng & Zhang (1999) present a model 
of multi-component X-ray emission from the rotation-powered pulsars, in which 
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hard thermal and soft thermal X-rays originate from the back-flow current of the 
outer gap. The directions of the soft X-rays are diffused, so that the resonant ICS 
may be effective due to proper incident angles. However, there is a criterion to 
determine whether a pulsar has an outer gap or not (Zhang & Cheng 1997, Cheng 
& Zhang 1999). The size of the outer gap, which is the ratio between the potential 
drop of the outer gap and the total potential of the open field lines, limited by the 
soft thermal X-rays from the pulsar surface, can be determined as, 

/ = 5.5P 26 / 21 Z?r 2 4/7 , (12) 

where P is the rotation period of pulsars in seconds, Bu the surface magnetic field 
strength in 10 12 G. Only when / < 1 can the outer gap exist. Unluckily, we find 
/ > 1 for all the seven pulsars, as presented in Table 2, and the outer gap may not 
exist. It seems that the criterion of Zhang & Cheng needs to be re-investigated, or 
other energy loss processes should be proposed. 

3). It should be emphasized that our result is model dependent. According 
to other radio emission models, e.g., relativistic plasma emission models (Melrose 
& Gedalin 1999 and references therein), some other interesting constrains on 
the relativistic particles or plasma may be obtained. But in this paper we just 
concentrate on the particle Lorentz factor so that only the Curvature radiation 
model and ICS model are employed. 

We are grateful to Dr. B. Zhang, Dr. J.L. Han and other members of the 
pulsar group for helpful discussion. This work is partly supported by NSF of 
China, the Climbing project, the National Key Basic Research Science Foundation 
of China, and the Research Fund for the Doctoral Program Higher Education. 
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Table 1: Parameters given by T91 and our results 



PSR B Freq. Range A A9 min n a ( r min 7CR,max 7ics,max r f^ 

' ' ' min ICR, max /ICS, max 





(GHz) 




(°) 




(°) 


(°) 


(km) 












0301+19 


0.1-2.7 


86 


0.9 


0.34 


31.9 


1.8 


59 


463 


3763 


4.4 


0.43 


0.092 


0329+54 


0.1-10.7 


1059 


19.8 


0.96 


30.8 


2.9 


169 


781 


2957 


2.2 


0.24 


0.068 


0525+21 


0.05-4.9 


90 


9.5 


0.47 


23.2 


0.7 


132 


762 


6105 


4.2 


0.28 


0.048 


1133+16 


0.026-10.7 


53 


4.4 


0.50 


51.3 


3.7 


135 


819 


4062 


2.8 


0.16 


0.031 


1237+25 


0.08-4.9 


79 


7.9 


0.52 


48.2 


0.9 


106 


621 


3419 


3.1 


0.31 


0.076 


2020+28 


0.1-14.8 


1103 


9.1 


1.08 


71.2 


3.6 


71 


666 


2374 


2.3 


0.22 


0.063 


2045-16 


0.08-4.9 


45 


7.9 


0.36 


36.7 


1.1 


132 


684 


4007 


2.7 


0.30 


0.078 



Table 2: numerical values of / 



PSR B 


P(s) 


B12 (G) 


/ 


0301+19 


1.388 


1.34 


6.98 


0329+54 


0.715 


1.21 


3.25 


0525+21 


3.746 


12.25 


6.74 


1133+16 


1.188 


2.11 


4.44 


1237+25 


1.382 


1.15 


7.58 


2020+28 


0.343 


0.81 


1.66 


2045-16 


1.962 


4.64 


5.27 
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Fig. 1. — The 7(7*) relationships for seven pulsars. The dot and solid curves 
represent the data calculated for curvature radiation and inverse Compton scattering 
respectively. 
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Fig. 2. — The effect of inner field lines to 7(7"). The curves marked by I are 
calculated for the last open field lines, and II for the inner field lines with A = 1.2. 
'CR' denotes curvature radiation and 'ICS' denotes inverse Compton scattering. 




Fig. 3. — The influence of inclination angle a and impact angle £ to j(r). The 
curves marked by I are calculated for (23.2°, 0.7°), II for (23.2°, 2.0°) and III (90°, 
0.7°). The marks of 'CR' and 'ICS' are the same as those in Fig.2. 



